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ABSTRACT 

We  derive  lower  bounds  for  the  norm  of  the  inverse  Vandermonde 
matrix  and  the  norm  of  certain  inverse  confluent  Vandermonde 
matrices.  They  supplement  upper  bounds  which  were  obtained  in 
previous  papers. 
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EXPLANATION 

The  sensitivity  of  a system  of  linear  alqebraic  equations  tcj 
small  perturbations  in  the  data  depends  in  large  measure  on  the 
magnitude  of  the  inverse  of  the  coefficient  matrix.  It  is  therefore 
of  interest  to  estimate  the  norm  (i.e.,  the  magnitude)  of  the 
inverse  of  a matrix.  We  do  this  here  for  the  Vandermonde  matrix 
and  certain  related  matrices,  which  occur  frequently  in  problems  of 
numerical  analysis,  providing  lower  bounds  for  the  norms  in  question. 
Upper  bounds,  and  exact  formulas  in  special  cases,  have  been  given 
previously . 
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ON  INVKRSKS  OF  VANDKKMONDE  ANU  CONFLUENT  VANDERMONUE  MATRICES  III 

Walter  Gautschi 


1.  Iiitrodiietion.  Norm  estimates  for  the  inverse  of  a Vandermonde  matrix,  or  the 
inverse  of  confluent  Vandermonde  matrices,  have  been  the  subject  of  several  previous 
pajjcrs  (1),  12],  (4).  The  emphasis  there  vas  on  upper  bounds  in  the  case  of  general 
complex  nodes,  or  identities  wlien  the  nodes  are  positive  [1],  12]  or  real  and  sj-mmetiic 
with  respect  to  the  origin  (4]  . VJe  now  wis)i  to  supplement  these  results  by  jiroviding 
lower  lx>unds  in  the  case  of  arbitrary  complex  nodes.  We  obtain  these  bounds  by  applying 
Jensen's  formula  in  the  tlioory  of  analytic  functions  to  appropriate  polynomials. 

2.  Jensen's  formula  for  polynoiiiials.  Given  a jxjlynomial 

(2.1)  p(z)=a  +asi"-+az'^,  a ^0, 

0 1 n n 

with  coinplex  coefficients  a , lot  5.  < . . • , denote  its  zeros  ordei  cd  ouch  that 

U i 2 n 

IcJ  1 IcJ  < 1 llj  1 1 < < |r,J  . 

Jensen's  f')imula,  apiilied  to  (2.1)  on  the  unit  circle,  then  gives  16} 
l^n"'r+l^r.2  t n | p (e"°)  |d0)  , 

I i 0 1 

licnee,  letting  M = max  |p(o  )|, 

0<0<2-'i 


a C C 
' n r+1  rn-2 


C I < M < I 
n - - . , 


L i h i I I maxd,  |c  |)  . 

li=0  " v=l 


Equality  in  (2.3)  liolds  if  and  only  if  a = a = ---=a  ,=0,  i.e.,  p(z)  - a z 

0 1 n~l  n 

Indeed,  if  p(z)  = then  (2.3)  (with  equality)  is  trivial.  Conversely,  if  we  have 
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equality  in  (2.3),  we  mur.t  have  equality  in  (2.2),  henre,  by  Jensen's  formula. 


iO 


p(o  ) I = M for  0 < 0 < 2t.  Since 


I , |2  r - i()<-f)0  V 

|l>(o  )|  - I aae  = I c.e 

)c,£=0  ^ X^.-n 


iA0 


is  a trigonometric  polynomial,  witli  coefficients 


'X  "li”);-X'  “-X 

K--«’ 


(the  coii\'(jntioii  — 0 if  M ^ 0 or  U ^ n is  used  here)  , it  can  be  constant  equal 


0 ] 


only  if 

C 

n 

‘'n-1  " 

• = = 0 and  Cp 

= m". 

The  first  condition. 

a a - 
n 0 

0, 

licnce  a,  - 
0 

0 (since  a ^ 0)  . 
n 

The 

second  condition. 

II 

1 

1 

0, 

then  gives 

a^^  - 0,  and  continuing 

in  this  manner,  we  find 

=;  . . . = 

a 

. = 0. 

n-1 


Invorsf  V.ujdonp^iidc  fr--^tvix.  V7e  denote  tlic  Vanderr.ionde  matrix  of  order  n by 

"i  1 ...  1 1 


(3.1) 


V (z)  - 
n 


"2 


n-1  n-1 

^ "-2 


n-1 


wh.  1 1 • /.  = 


^1 ' ^'2' ’ ‘ ' ^'n'  in  a vector  of  n complex  numbers,  calleS  "nodes".  If  tl>e 


nodes  are  mutually  di.stim  t,  then  has  an  inverse,  v;)iicli  wo  denote  by 


(3.2) 


A , )i-l 


Wo  air  interested  in  t.lio  C -norm  of  ('H.2), 


V ^ 11^  = I 


Xu'  • 


l_^A<n  u-A 

Th.-ort^m  ,1.1.  If  distinct  complex  numbers , and  n > 1,  th 


(3.^) 


V (z) 
' n 


n max (1 , z ) 
max  I ' 


l<X<n  v-1  z. 


- 2 


Proof.  We  recall  lH  that  the  elements  u in  (3.31  arc  the  coefficients  of  the 

Aji 

fundamental  Lagranqc  intcriK>lation  polynomials  associated  with  the  nodes  z , 

V 


(3.4) 


n 

v=l  “X 


= u + u z + 

z.  - z XI  X2 


+ u.  z 
Xn 


•n-1 


Applying  (2.3)  and  t)ie  remark  following  (2.3)  to  the  polynomial  of  degree  n-1  in  (3.4), 
we  f i nd 


(3.5) 


^1  ^ y T^: 

U=1  v^X  ' X 


1 


7~T  maxd,  |z  I ) . 
v^X 


It  ^0  index  X for  wliicli  the  right-hand  expression  in  (3.5)  attains  its 

n 

maximum,  then  that  maximum  is  less  than  ^ |u  |,  •fionce  less  or  equal  than 
n U=1 

max  '[  |u  |.  This  cstat>lishes  (3.3)  and  ]irovcs  Theorem  3.1. 
l<X<n  „=1 

The  lower  bound  in  (3.3)  supplfii\ent.';  tlic  upper  (attainable)  bound  in  [1],  which  is 
of  the  sai;ic  form  as  (3.3)  except  tliat  the  f^-norm  of  the  2-vectors  |l,z^)  in  tlic  numerator 
factors  is  leplaceJ  Iry  tlie  £^-norm. 


4.  Inverse  conf lucid  Vaiidermor.dn  matrices.  TliO  technique  used  in  tlic  proof  of 
Theorem  3.1  can  bo  adaptc^d  to  confluent  Vandermonde  matrices.  We  illustrate  this  with 
the  particular  matrix 


(4.1)  lh_  (z)  = 
2n 


2n-l  2n-l 

Zl  z^ 

considered  previously  in  11),  12). 


0 

1 

?z. 


2z. 


2z 


2n-l  2n-2  2n-2  ,,  2n-2 

z (2n-l)z,  (2n-l)z.  ...  (2n-l)z 

n 1 2 n 


'^oorem  4.1.  If  Zj,z^,...,z^  ajre  mutual  1y  distinct  complex  numbers , and  n > 1,  then 

2 


(4.2) 


cl  max(l,  Ik  |)  \ 

|U,;(Z)  IL  > max  b^rj  -K  --tV 
lf^X5_n  v=l  \ ' X v'  / 

v/X  ' ' 


-3- 


whc're  ^ the  lai <|<T  of  t In^  two  quant itios 

(4.3)  b*^^  = max(l  , |zj^  I ) , b^^^  = max|2j  J,  !/(?• 

Proof . We  have  12] 


U 


-1 

2n 


[V,  1, 
Xp 


w = |w  ) , 
Xp 


wl)ere 


2n 


tj(z)|l  - 2fj^(Zj^)(z  - Zj))  = I z*'  ^ 


(4.4) 


p=.l 


2n 


Xp 


1 < X < n , 


f^(z)(z  - z^)  = I w z 
p=l 


P-1 


and  donotc.s  the  fundamontal  Lagrange  interpolation  polynomial  in  (3.4).  Applying 

(2.3)  to  tlic  polynomials  in  (4.4),  and  talking  note  of  tlio  rem.u)'.  following  (2.3),  one  finds 


where  delincal  in  (4.3). 

Ucnot  iiiq 

the  products 

] \ on  the  right  by 

It,,  and  observing  that  III'  ^ II  = max  (max 
X :■  M 2n  “■  ^ 

2n 

2n 

X 

1 

P=1 

l'’xpl' 

max  y |w, 

1 ' Xp 

X |i=l 

|),  an  argument  similar 

to  tlio  one  after  (3.5)  will  show  that 

X <1: 

1%^.  11^ 

1 (2)  II 

0'  ^ "x  ■ II 

U ^ II  for  all 

X = 1,2, ...,n,  hf  nee  max  Or  ,b^  ”x  ^ ' 

1 "2; 

II  _ for  all  X 1 . 

, 2,...,n.  This  proves 

Theorem  4.1. 

The  lower  l>onnd  in  (4.2)  r.upj'li. nont s the  (attainable)  upi'f-r  bound  in  (2),  whicli  is 
of  tlie  same  fo'm  as  (4.2)  I'Xo.pt  Diat  tlii>  t^-norm  of  the  2-vectors  (],z^)  in  the  numei  ator 
factors,  and  tlie  t^-norms  defining  ptifl  replaced  by  the  respiective 

l^-noims.  Tn  the  case  of  positive  nodes,  another  (usually  sharper)  lower  bound  can 


bo  found  in  13,  Theonm  2.1). 


5 . Examplor. . 

Cxam)'lo  5.1  (roots  of  unity),  v = l,2,...,n. 

In  view  of 

j n ,U-1 

" n j^l  ' ^ ' 

VO  obtain  fron  (3.4),  and  from  (4.4)  after  a little  computation, 

(5.1)  ||v"^z)|l  =1,  ||u'^(z)|l  =2--. 

n " 2n  ' " n 

Tile  lower  bounds  in  (3.3)  and  (4.2)  both  evaluate  to  1/n,  while  tlie  uppei'  bounds  in  (1)  , 
12)  are  2’^  Vn  and  (2n  - 1)4^  V»^,  respectively. 

Example  5.2  (rool:s  of  unity  on  half-circle).  z^  = (v-l)/N^  ^ _ y^2,...,n,  wlicre 

n = (N/2)  + 1, 

The  true  norms  of  V ^ and  U as  well  as  the  lower  bounds  of  Theorems  3.1  and 
n 2n 

(2) 

4.1  and  Ujc  upper  bounds  in  (1),  12)  are  shown  in  Taiile  5.1  for  N = 5(5)20  . It  is 

interestincj  to  note  how  deletion  of  the  roots  of  unity  on  a half-circle  results  in 

sulis'.  inticil  ly  lar'jc)  v.aluos  of  llv  ^11  and  lllJ  Ml 

n n "■»  " 2n  »■ 


ii"2;ii» 

N 

n 

lov..-r 

true 

ui>j  er 

lower 

true 

upper 

5 

3 

7.24 (-1) 

1.89 

2.89 

1.57 

1.79(1 ) 

4.19(1) 

10 

1.17 

1 .47(1) 

3.75(1) 

8.29 

2.36(3) 

1 .56(4) 

15 

8 

4.25 

2.03(2) 

5.45(2) 

1.46(2) 

6.18(5) 

4.48  (6) 

20 

1 1 

1.17(1) 

2.70(3) 

1 .20(4) 

1.52(3) 

1 .59(8) 

3.03(9) 

Table  5.1.  Norm  estimates  for  Example  5.2. 

2 n 

Examide_5^3.  e (z  ) = 0,  v = l,2,...,n,  where  e (z)  = 1 + z + --+•••+  — . 
~ V n 2 ! n I 

Ilsiivj  (he  zeros  of  e^,  tal'iilated  in  15),  we  obtain  the  results  in  Talile  5.2. 

(27 


The  inlefjers  in  j'arcntb.eses  indicate  exiMnents  of  10. 


1.12 

2.00 

3.93 

4.76 

2.21  (1) 

7.90(1) 

2.44 

5.22 

1.45(1) 

4.45(1) 

2.52(2) 

1.96(3) 

7.45 

1.69(1) 

5.71 (1) 

6.08(2) 

3.69(3) 

4.22(4) 

2.27(1) 

5.36(1) 

2.02(2) 

7.54 (3) 

4.79(4) 

6.84 (5) 

Tdbli  5.2.  Nyim  estinatos  for  Kxa.T.plo  5.3. 

Example  5.1.  c (z  ) ■■■  0,  Im  z > 0,  u = l,2,...,n,  where-  n = . 

— N V V — L 2 J 


N 

n 

lower 

true 

upper 

5 

1.62 

2.74 

3.13 

10 

■■ 

5.71 

1.09(1) 

1 .27  (1 ) 

15 

8 

3 . 07  ( 1 ) 

6.82(1) 

8.63(1) 

20 

10 

1.60(2) 

3.60(2) 

4.49(2) 

6.69 
1.3H(2) 
5.71 (3) 
l.H8(5) 


2.51 (1) 
6.15(2) 
3.24 (4) 
1.07(6) 


upper 

3.29(1) 

8.35(2) 

5.19(4) 

1.66(6) 


Table'  5.3.  Nr.rin  eolinati-s  f ni  r.xair.pli.'  5.4. 

-Similarly  .ui  i ti  I'.x.iinpli-  5.2,  th  li'lion  of  Uk-  zeros  in  Uie  )ov;or  lialf-plaiie  has  tlic 
effc.-t  of  i n -I  I -i:;i  ii'i  I lio  noin-i  of  V ^ and  I',.'. 
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